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We have made an angle-resolved photoemission study of a one-dimensional (ID) Mott-Hubbard 
insulator Nao.geV^Os and found that the spectra of the V 3d lower Hubbard band are strongly 
dependent on the temperature. We have calculated the one-particle spectral function of the one- 
dimensional t-J model at finite temperatures by exact diagonalization and compared them with 
the experimental results. Good overall agreement is obtained between experiment and theory. The 
strong finite temperature effects are discussed in terms of the existence of the "Fermi surface" of 
the spinon band. 

PACS numbers: 71.27.+a, 79.60. Bm, 71.10.Pm, 75.20.Ck 
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The most striking and non-trivial theoretical predic- 
tion for ID strongly correlated systems is the spin-charge 
separation the degrees of freedom of an electron 
are decoupled into elementary excitations of spin and 
charge called "spinon" and "holon" , respectively. Angle- 
resolved photoemission spectroscopy (ARPES) is a pow- 
erful technique to study the spin-charge separation: Kim 
et al. Q performed a pioneering ARPES work on the 
ID charge-transfer insulator SrCuC>2 and found that the 
spectra agree well with the theoretical one-particle spec- 
tra of the ID t-J model with realistic parameters, iden- 
tifying such spinon and holon excitations. We subse- 
quently investigated whether the same scenario is appli- 
cable for a ID Mott-Hubbard type insulator NaV^Os ||. 

Strongly-correlated systems are often characterized by 
the presence of a characteristic low energy scale in spite 
of the large energy scales of bare interaction strengths. 
For the Hubbard model, the relevant low energy scale is 
set by the superexchange interaction J ~ it 2 /U (<C t, U) 
rather than the bare interaction parameters, the transfer 
integral t and the on-site Coulomb repulsion U. Since the 
photoemission spectrum is a projection of the initial state 
onto the set of final states, drastic finite-temperature ef- 
fects may be expected for a temperature change of the 
order of such a characteristic low energy scale. Finite- 
temperature effects can be particularly drastic in ID sys- 
tems j| because of the existence of the "Fermi surface" of 
spinon excitations. In this Letter, we present the result 
of a temperature-dependent ARPES study on NaV^Cv 
While severe charging effects had previously prevented 
the measurements of NaV^Os below ~ 300 K Q, more 
conductive Nao.geV^Os || enabled us to obtain ARPES 
spectra at temperatures as low as ~ 120 K. Dispersive 



features of V 3 c? character in the lower Hubbard band 
were found to be dramatically dependent on the temper- 
ature. In addition, we have made a comparison with the- 
ory and confirmed that the observed finite-temperature 
effect is due to strong correlation effect rather than the 
simple thermal broadening. 

Recently NaV^Oa has widely aroused much interest as 
a quasi-lD quantum-spin system. Crystalographically, it 
is still controversial whether the V atoms in VO5 pyra- 
mids running along the &-axis of NaV^Os is mixed- valent 
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In the former case, this compound can be no doubt re- 
garded as a half- filled chain , while in the latter case it 
is viewed as a quarter-filled ladder system |^|,||. At low 
temperatures around or below its spin-Pcicrls-like (SP) 
transition temperature T$p ~ 34 K 0, the difference 
between those two models is significant in terms of the 
charge ordering pattern. On the other hand, the mag- 
netic susceptibility x(T) well above T$p is successfully 
fitted to the Bonner-Fischer curve with J ~ 560 K 0, 
indicating that NaV^Os behaves as a good ID antiferro- 
magnetic Heisenberg chain in this temperature region. In 
fact, it is theoretically supported that it can be mapped 
on to the ID Heisenberg chain ||. Except that the SP 
transition is suppressed ||, the Na-deficient Nao.96V205 
has almost the same magnetic properties as NaV^Os [ (To] . 
Though remaining an insulator, Nao.geV^Os is more con- 
ductive than NaV^Os due to the doped holes. 

Single crystals of Nao.geV^C^ were prepared as re- 
ported in Ref. They could be easily cleaved par- 
allel to the ab plane. The ARPES measurements were 
made using the He I resonance line (hv = 21.2 eV) and a 
hemi-spherical analyzer with angular resolution ±1° and 
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energy resolution 80-100 meV. The measurement tem- 
perature ranged from T = 120 K (= 0.21J) to room 
temperature 300 K (= 0.54 J). The measurements were 
performed for several in situ cleaves, for which we care- 
fully cycled the temperatures of cleavage and measure- 
ments in order to exclude any extrinsic effects such as 
surface degradation and contamination. 

Before discussing the V 3d band features of our main 
interest, we would like to mention that the O 2p band 
structure was found very anisotropic: ARPES spectra 
with momentum parallel to the &-axis (fc || 6) show rich 
dispersive features while those with momentum parallel 
to the a-axis (k _L b) have no dispersions, supporting the 
one-dimensionality of this compound. These results well 
agree with those of stoichiometric NaV2C>5 j|, indicating 
that the Na deficiency has no appreciable influence on 
the O 2p band structure. Furthermore, the O 2p spectra 
have no obvious temperature-dependent changes, making 
a clear contrast to the remarkable changes of the V 3d 
band described below. 

Figures [j] (a) and (b) show the geometry of the mea- 
surements. By scanning the momentum in the Brillouin 
zone (BZ) as indicated, we obtained ARPES spectra 
along the k _L b and "k || b" cuts as shown in Figs. |l| 
(c) and (d), respectively. Hereafter, k a and kb are the 
momenta along the a- and b- axes in units of the in- 
verse of the corresponding lattice constant, respectively. 
All the spectra have been normalized to the area of the 
lower Hubbard band of V 3d origin between the binding 
energies (Eb) of 0.0 and 2.5 eV after subtracting the inte- 
gral background. We notice that the k Lb cut [Fig. |l|(c)] 
shows a clear temperature dependence and a weak angle 
dependence except for around the T point (6 — <p = 0°). 
In addition, the absolute intensity at the T point was 
found to be weaker than half of those at <f> > 10°. These 
observations can be explained by the facts that the occu- 
pied V 3d orbital has xy symmetry lying approximately 
in the ofe-plane [jl2|] and therefore that normal emission 
from this orbital is forbidden due to selection rules 
We therefore conclude that the temperature dependence 
results from the intrinsic finite temperature effects of the 
V 3d xy band |]l4|| and that the momentum dependence 
along the k _L b direction is due to the matrix element 
effects of the d xy orbital. 

Figure [l] (d) shows the result for the "fc || 6" cut. In or- 
der to avoid the T point where the matrix- element effects 
prohibit emission from d xy , the k a value was slightly off- 
set from the 6-axis as shown in Fig. [I] (b) . The obtained 
120 K spectra show rich dispersing features. In going 
from kb = to 7r/2, a single peak centered at Eb = 0.9 
eV is split into two features: the splitting becomes largest 
at kb = ir/2 with the two features located at Eb ~ 0.7 
and ~ 1.4 eV. The 0.7 eV peak then decreases in intensity 
in going from kb = tt/2 to n and only a single broad peak 
is left at Eb ~ 1.1 eV at the BZ boundary kb ~ tt. The 
fch-dependence of the spectra between kb = and 2ir is 
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FIG. 1. (a) Definition of the take-off angles cj> and 9 of 
photoelectrons. (b) Measured cuts in the the Brillouin zone, 
(c) Spectra for the k _!_ b cut and (d) those for the "k |j 6" cut. 
Solid and dashed curves show spectra taken at 120 K and 300 
K, respectively. Note that each spectrum is normalized to 
its area and that the absolute intensities of the spectra near 
6 = <j> = 0° are much weaker than the others. 

almost symmetric with respect to kb — tt, which excludes 
significant changes in the photoemission matrix elements 
between the first and second BZ's. By contrast, the 300 K 
spectra, which well agree with the previous report ||,[l5| , 
show less pronounced features than the 120 K spectra. 
The spectra at kb = and 7r become a broader peak with 
a longer tail towards high binding energies. As for the 
spectra at kb ~ 0.57T, the peak located at Eb ~ 0.7 eV 
becomes weaker and that at ~ 1.4 eV stronger in going 
from 120 K to 300 K. We also performed ARPES mea- 
surements at 200 K and confirmed that the changes are 
gradual as a function of temperature. These observa- 
tions are more clearly recognizable in the intensity plots 
(b) and (c) of Fig. §. 

By noting that the finite-temperature effects strongly 
depend on the momentum (k = kb), they are clearly 
not due to simple thermal broadening or charging ef- 
fects. The following two points may be remarked: (i) 
in the temperature range studied here, there is no phase 
transition in this compound that may give rise to such 
a dramatic change; (ii) the energy scale of the spectral 
change is not of order kT ~ 0.03 eV but of order ~ 1 eV. 
These phenomena are obviously beyond the conventional 
band picture and would reflect strong correlation effects. 

In order to interpret the above observations quantita- 
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FIG. 2. Comparison between the experimental and the- 
oretical spectra. Panel (a) shows the ARPES spectra of 
Nao.96V20s with the intensity plots in panels (b) and (c). 
Panels (d), (e) and (f) show the theoretical counterparts cal- 
culated for the half-filled 14-site t-J model. Here J = 0.05 
eV. 



tively, we have adopted the ID t-J model to calculate 
the one-particle spectral function A(k,u) at finite tem- 
peratures by the exact diagonalization method Jlq| . As 
the behavior of x(T) well above T$p indicates, the t-J 
model is valid as one of the simplest models to describe 
this system in the temperature region considered. It may 
be applicable not only to the half-filled chain case 
but also to the quarter-filled ladder case where each d 
electron is localized in a rung of two V atoms In 
the ladder case, only the half-filled bonding band is taken 
into account, because no electrons are in the antibonding 
band [||. The quarter-filled ladder can be mapped on to 
the half-filled chain by taking only the half-filled bond- 
ing band with the empty antibonding band neglected and 
by an appropriate modification of the parameters t and 
J HJD]]. I n addition, the t-J model has the advantage 
over a more realistic Hubbard model in that the former 
can treat larger clusters, which is crucial for the discus- 
sion of finite temperature effects. We have calculated 
A(k,u); 
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where Z = Y2i e l3Ei is the partition function. Results 
for a half-filled 14-site t-J cluster with J It = 1/3 at 
T = 0, J/4 and J/2 are shown in Fig. | §J. A(k,uj) at 
T = can be intuitively interpreted as a convolution of 
spinon and holon excitations, whose dispersional widths 
are ~ J and ~ 2t, respectively. In this scenario, in the 
ground state, the holon band is empty while the spinon 
band is half-filled up to the Fermi momentum k = w/2 |Q. 
The lineshape of the "spinon branch" (see Fig. ||) is de- 
termined by the band edge singularity of the holon band 
whereas that of the "holon branch" is determined by the 
Fermi edge singularity of the spinon band Jl9) . At finite 
temperatures, the spectral weight of the spinon branch 
at u/t > 1 (u>/t < —2) is transferred from < k < tt/2 
(tt/2 < k < 7r) to tt/2 <fc<7r(0<fc< tt/2). At the 
same time, the intensity of the holon branch decreases. 
In fact, spectral weight is transferred from the spinon 
branch to a wide energy region, making the spectral fea- 
tures less pronounced and more symmetric with respect 
to k = n/2. In fact, the singularity of the holon branch 
due to the existence of the spinon Fermi surface is easily 
smeared out over the entire energy range of 4t (3> T ~ J) 
at finite temperatures of order J 01. 
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FIG. 3. Spectral function of the half-filled 14-site t-J model 
with J/t = 1/3 at various temperatures. The thick solid 
and dashed curves indicate the singularities of the spinon and 
holon branches, respectively. 

In Fig. U we show a comparison between the experi- 
mental and theoretical spectra. Figure || (a) shows exper- 
imental spectra for selected momentum values and Figs. ^ 
(b) and (c) show the intensity plot in the E-k plane on the 
grey scale. Correspondingly, we show in Figs. ^ (d)-(f) 
the theoretical spectra at T — J/4 ~ 150 K (solid lines) 
and T = J/2 ~ 300 K (dashed lines) after Gaussian 
broadening with the width 2.5t and their intensity plots. 
Here, we have assumed that t = 3J = 0.15 eV, which 
is plausible because J = U 2 /U ~ 0.05 eV and U = 2-4 
eV in typical vanadium oxides pfj. Comparing (a)-(c) 



3 



and (d)-(f), overall agreement is satisfactory. In the low 
temperature region (T — 120 K or J/4), the experimen- 
tal shift of the peak position between kb = 0.0 and 0.957T 
may be attributed to the existence of the spinon branch, 
resulting in the asymmetry of the spectra with respect to 
k = 7r/2. Besides, between ki — 0.32-7T and 0.727T, there 
are two dispersing features which may be assigned to the 
two holon branches as reproduced in the theoretical spec- 
tra. These findings are also substantiated in the compar- 
ison between (b) and (e). Unlike the ID cuprates, where 
the intense O 2p structure obscures the higher binding 
energy part of the holon and spinon branches [0, the 
whole structure of theoretical A(k, oS) can be compared 
with the experimental results of this compound. 

In the high temperature region (T = 300 K or J/2), 
both results become broader and more symmetric with 
respect to k = ir/2. As a result, the tendency of the 
experimental spectral weight redistribution due to finite 
temperature is grossly reproduced by the theory as seen 
in the intensity plots in Fig. ||. Around k ~ 0, agreement 
is quite excellent. To this extent, the experimental finite 
temperature effects may be attributed to the existence of 
the spinon Fermi surface, which theoretically causes the 
dramatic spectral redistribution over the entire E-k space 
with changing temperature. To be more precise, however, 
there exist some discrepancies between theory and exper- 
iment. When the temperature is increased from 120 K 
to 300 K, the spectra change more dramatically than the 
theoretical prediction. Furthermore, while the temper- 
ature dependence is rather well simulated by theory at 
feft ~ 0-0.32-7T, at fcb ~ 0.57T the feature around Eb = 0.7 
eV in experiment loses much of its spectral weight in go- 
ing from 120 K to 300 K, in disagreement with theory. In 
addition, around fc& ~ ir in going from 120 K to 300 K, 
the experimental spectra lose spectral weight at Eb < 0.8 
eV and a longer tail develops on the high binding energy 
side unlike the theoretical spectra. At low temperatures, 
where the decay of a photohole is probably dominated by 
purely electronic mechanisms while other decay channels 
may become available at higher temperatures. As a can- 
didate, we may consider electron-phonon interaction and 
possible charge disorder, which may become important 
at higher temperatures. The difference between the t-J 
model and the Hubbard model as well as the degener- 
acy of the V 3d orbitals might be another origin for the 
discrepancy. 

In conclusion, we have made an ARPES study of 
Nao.96V205 by changing the temperature and found that 
a strong spectral weight redistribution occurs in the lower 
Hubbard band. Also we have calculated the one-particle 
spectral function of the ID t-J model at finite tempera- 
tures by the exact diagonalization method. The overall 
agreement between theory and experiment implies that 
the spin-charge separation picture is valid in this system. 
Although they are more drastic than the theoretical pre- 
diction, the experimental finite temperature effects have 



been partly explained by the theory, which may be ex- 
pressed as the "Fermi surface" effect of the spinon band. 
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